Isoscaling in the Lattice Gas Model 
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The isoscaling behavior is investigated using the isotopic/isobaric yields from the equilibrated 
thermal source which is prepared by the lattice gas model (LGM) for lighter systems with A = 36. 
The isoscaling parameters a and -j3 are observed to drop with temperature for the LGM with the 
asymmetric nucleon-nucleon potential. However, the isoscaling parameters do not show temperature 
dependence for the LGM with the symmetric nucleon-nucleon potential. The relative neutron or 
proton density shows a nearly linear relation with the N/Z ( neutron to proton ratio ) of system. 

PACS numbers: 25.70.Pq, 24.10.Pa, 05.70. Jk 



Isoscaling has been observed in a variety of reactions 
under the conditions of statistical emission and equal 
temperature recently by Tsang et al P, 0, ■ This kind 
of scaling means that the ratio R,2i(N,Z) of the yields of 
a given fragment (N,Z) exhibits an exponential depen- 
dence on N and Z when these fragments are produced 
in two reactions with different isospin asymmetry, but at 
the same temperature. Experimentally the isoscaling has 
been explored in various reaction mechanisms, ranging 
from the evaporation [| , fission 0, [j| and deep inelastic 
reaction at low energies to the projectile fragmentation 
la, H and multi-fragmentation at intermediate energy 
pUsH- While, the isoscaling has been extensively ex- 
amined in different theoretical frameworks, ranging from 
dynamical model, such as Anti-symmetrical Molecular 
Dynamics model and BUU model 0, to statistical 
models, such as Expansion Emission Source Model and 
statistical multi- fragmentation model |2ll3l lllLll2| . From 
all these reaction mechanisms and models, it looks that 
isoscaling is a robust probe to relate with the symmetri- 
cal term of the nuclear equation of state. 

Typically, the investigations of isoscaling focused on 
yields of light fragments with Z—2-8 originating from de- 
excitation of massive hot systems produced using reac- 
tions of mass symmetric projectile and target at inter- 
mediate energies, such as 112,124 <SVi + U2 ' 12i Sn in Michi- 
gan State University (MSU) data 0, 0, Q or by reac- 
tions of hig h- energy light particle with massive target 
nucleus . In a recent article @ , the isoscaling us- 

ing the heavy projectile residue from the reactions of 25 
MeV/nucleon 86 Kr projectiles with 124 Sn, 112 Sn and 64 Ni, 
58 Ni targets which was performed at Texas A&M Univer- 
sity (TAMU) and the isoscaling phenomenon on the full 
sample of fragments emitted by the hot thermally equi- 
librated quasi-projectiles with mass A = 20-30 are also 
reported 0. 

In this study, we present an isoscaling analysis for the 
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light fragments from thermal sources which are produced 
by the lattice gas model. Instead of the fixed charge 
number of the reaction system in MSU data or TAMU 
data, here we fixed the source mass for lighter system, 
but changing the charge number and neutron number. 
The isospin fractionation is also observed via the relative 
free neutron and free proton density which is obtained 
by the isoscaling parameters. 

A thermally equilibrated system undergoing statistical 
decay can be, within grand-canonical approach, charac- 
terized by a yield of fragments with neutron and proton 
numbers N and Z 

Hill: 

Y(iV,Z) = F(Y,Z)exp^exp(^i + ^) (1) 

where F(N, Z) represents contribution due to the sec- 
ondary decay from particle stable and unstable states to 
the ground state; /i n and \x v are the free neutron and 
proton chemical potentials; B(N, Z) is the ground state 
binding energy of the fragment, and T is the tempera- 
ture. 

The ratio of the isotope yields from two different 
systems, having similar excitation energies and similar 
masses, but differing only in N/Z, cancels out the effect of 
secondary decay and provides information about the ex- 
cited primary fragments 0. Within the grand-canonical 
approximation ( Eq.lJTJ j, the ratio Y 2 {N,Z)/Y l (N,Z) 
assumes the form 

R 2 i(N,Z) = Y 2 (N, Z)/Y 1 (N,Z) = Cexp(aN + pZ) 

(2) 

with a = A/i„/T and (3 = A/i p /T, with A/i„ and A/i p 
being the differences in the free neutron and proton chem- 
ical potentials of the fragmenting systems. C is an overall 
normalization constant. 

The tool we will use here is the isospin dependent 
lattice gas model (LGM). The lattice gas model was 
developed to describe the liquid-gas phase transition 
for atomic system by Lee and Yang [0|. The same 
model has already been applied to nuclear physics for 
isospin symmetrical systems in the grand-canonical en- 
semble |l7| with a sampling of the canonical ensemble 
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ei Ens mm mm, and also for isospin asymmet- 
rical nuclear matter in the mean field approximation [25| . 
Here we will make a brief description for the models. 

In the lattice gas model, A (= N + Z) nucleons with 
an occupation number s which is defined s = 1 (-1) for 
a proton (neutron) or s — for a vacancy, are placed on 
the L sites of lattice. Nucleons in the nearest neighboring 
sites interact with an energy e SiS .. The hamiltonian is 
written as 



A P 2 

i— 1 i<j 



(3) 



In order to investigate the symmetrical term of nuclear 
potential in this model, we use two sets of parameters: 
one is an attractive potential constant e SiSj between the 
neutron and protons but no interaction between like nu- 
cleon, i.e. proton and proton or neutron and neutron, 
namely 



O.MeV, 



-5.33MeF. 



(4) 



This potential results in an asymmetrical potential 
among different kind of nucleons, hence it is an isospin 
dependent potential. For simplicity we call the calcula- 
tion with this potential as isoLGM thereafter. Another 
set is the same interaction constant between like nucleons 
or unlike nucleons, i.e. 



-hMMeV. 



(5) 



In this case, the nucleon potential is symmetrical among 
all nucleons, i.e. isospin independent potential. For sim- 
plicity, we call the calculation with Eq.(5) as noisoLGM 
thereafter. In this work, mostly we use isoLGM to ex- 
plore isoscaling behavior, but we will also use noisoLGM 
to compare the isoscaling results. 

In the LGM simulation, a thrcc-dimension cubic lattice 
with L sites is used. The freeze-out density of disassem- 
bling system is assumed to be pf = j;P0i where po is the 
normal nuclear density. The disassembly of the system is 
to be calculated at pf , beyond which nucleons are too far 
apart to interact. Nucleons are put into lattice by Monte 
Carlo Metropolis sampling. Once the nucleons have been 
placed we also ascribe to each of them a momentum by 
Monte Carlo samplings of Maxwell-Boltzmann distribu- 
tion. 

Once this is done the LGM immediately gives the clus- 
ter distribution using the rule that two nucleons are part 
of the same cluster if P^/2p, — e SiS SiSj < 0. This method 
is similar to the Coniglio-Klein's prescription 26] in con- 
densed matter physics. 

In this paper we choose the small size nuclei with A 
= 36 as emission sources. Four isotonic sources, namely 
36 C*a, 36 Ar, 36 S and 36 Si, corresponding to N/Z = 0.8, 
1.0, 1.25 and 1.57, respectively, are simulated. In all 
cases, the freeze-out density pf is chosen to be 0.563 po, 
which corresponds to 4 3 cubic lattice is used. 10000 
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FIG. 1: (Color online) Isoscaling behavior of R.2i(N) (the upper 
panel) and R2i(Z) (the lower panel) at T = 5.0 MeV. 
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FIG. 2: (Color online) Isoscaling parameter a and /3 extracted in 
fixed Z from Fig. ^a) and N from Fig. nTb - ). respectively. 



events were simulated for each T which ensures good 
statistics for results. 

As an example, in Fig. ^ we present the isotopic scal- 
ing (the upper panel) and the isobaric scaling (the lower 
panel) from the isospin dependent lattice gas model for 
36 S to 36 Ca, respectively, at T = 5.0 McV. The charged 
particles with Z<8 and 2 < N < 9 have been accumulated 
to perform the ratios. There exist good linear behavior 
in semi- log plot. In order to extract the isoscaling param- 
eters a and /3, we use R 2 i{N) = Cexp(aN) to obtain a 
for a given Z (Fig. ^a)) and use R.2i(Z) — Cexp(f3Z) to 
obtain /3 for a given N (Fig. ^b)), respectively. Fig. [3 
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FIG. 3: (Color online) The absolute value of the isoscaling param- 
eters a and f} as a function of T from the yield ratios from the 
sources of 36 S and 36 Ca. The scattering points are the calcula- 
tion results for LGM with asymmetrical nucleon-nucleon potential 
(Eq.(4)) and the lines represent the results for the LGM with the 
same nucleon-nucleon potential for unlike and like nucleon pairs 
(Eq.(5)). See text for details. 



shows the extracted a or (3 versus Z or N . Obviously, a 
or (3 keeps the same value in the wide Z or N range and 
their absolute values are almost the same which is due 
to the absence of Coulomb interaction in the lattice gas 
model. 

Fig. shows the temperature dependence of the abso- 
lute values of isoscaling parameters a and (3. Here the 
scattering points are the calculation results for LGM with 
asymmetrical nucleon-nucleon potential and the lines 
represent the results for the LGM with the same nucleon- 
nucleon potential for unlike and like nucleon pairs. A 
decreasing trend of the values is clearly seen when the 
asymmetrical nucleon-nucleon potential is used in LGM, 
which indicates that the isospin dependence of the frag- 
ment yields becomes weak with increasing temperature. 
Considering that a = Afi n /T and (3 ~ A/ip/T, we can 
deduce the differences in free neutron and proton chemi- 
cal potentials of the fragmenting systems which is shown 
in Fig. 2Ja). Within the error bars, it looks like that 
Afi n and A/i p keep constant in the case of asymmetri- 
cal potential (Eq.(4)). However, a slight and wide valley 
can be also identified around T = 5 MeV as well as a 
slight kink shows in SMM model for Sn systems [Tlj . 
For our system, this valley may be related to the liquid 
gas phase change for the similar system around 5 MeV 
in the same model calculation as well as the data p7j . 
This valley becomes more obvious if we plot the values 
of (A/i„ — Afj, p )j2 as a function of temperature as sug- 
gested in Rcf. 0, E3. Fig. gfb) shows that a turning 
point seems to occur around 5 MeV. Of course, the error 
bars look larger for lower and higher temperatures, be- 
cause of the poor statistics for diverse cluster species due 
to the evaporation mechanism in low T and the vapor- 
ization mechanism in higher T. Around this point, an 
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FIG. 4: (a) (Color online) The absolute value of the difference 
of neutron and proton chemical potential (Afi n and |A/i p |) as a 
function of T. (b) (Afi n — Afi p )/2 as a function of temperature. 
Solid points are calculated results and line is its second polynomial 
fit. See text for details. 



apparent critical behavior has been observed in the dis- 
assembly of hot nuclei with A~36 in TAMU-NIMROD 
experimental data and model calculations :2?l by a wide 
variety of observables: such as the maximal fluctuations, 
critical exponent analysis and fragment topological struc- 
ture, namely nuclear Zipf law of Ma |2g and the heaviest 
and second heaviest correlation etc. Hence, this turning 
point of the difference of neutron and proton chemical 
potential might give an additional evidence of the chem- 
ical phase separation when the liquid gas phase transi- 
tion occurs. Experimentally, this kind of turning point 
has been observed recently for the quasi-projectile from 
the peripheral collisions of 28 Si + 112,124 SVi at 30 and 50 
MeV/nucleon in another TAMU data and can be under- 
stood as a signal of the onset of separation into isospin 
asymmetric dilute and isospin symmetric dense phase in 
a recent paper by Veselsky et al 0. 

Further, we shall investigate the effect of the asymmet- 
rical nucleon-nucleon potential in the lattice gas model 
on the isoscalig behavior. We take the same potential 
between like nucleons and unlike nucleons in the lattice 
gas model (Eq.(5)) to compare with LGM with isospin 
asymmetrical nucleon-nucleon potential (Eq.(4)). In this 
case, we observe that temperature dependence of a and 



4 




Finally, we test the relative free neutron density and 
proton density, which can be deduced by the equations: 



p n = exp(a'), 
p p = exp(P'). 



(6) 
(7) 
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FIG. 5: (Color online) (a) The relative free neutron (solid symbols) 
or proton (open symbols) density as a function of the N/Z of sys- 
tem at T = 5 MeV for the isoLGM case (squares) and noisoLGM 
case (circles), (b) The solid and dashed line represents the ra- 
tio Pn/Pp for isoLGM case and noisoLGM case, respectively. The 
dotted line shows the initial pn/Pp assuming neutrons and protons 
homogeneously distributed in a volume proportional to the nucleon 
number. 



(3 vanishes and their absolute values are smaller than the 
isoLGM case as shown by the lines in Fig. The insensi- 
tivity of a and (3 to temperature means that the isoscal- 
ing of fragment yields only stems from the initial N/Z 
of the system. This is possible since the nucleon-nucleon 
potential is the same for like nucleons and unlike nuclc- 
ons in this noisoLGM case. When temperature increases, 
the yield of light clusters also increase, but their isospin 
population does not change between the different cluster 
species, which results that no dependence for the isotopic 
or isobaric ratio on temperature. In this case, the system 
looks like a single component fluid, the chemical phase 
separation can not occur. However, for the isoLGM, the 
nucleon-nucleon potential is asymmetrical between like 
pair and unlike pair, it will give an additional isoscaling 
contribution due to the isospin dependent equation of 
state except from the trivial contribution from the initial 
difference of the isospin of the source, which results in 
a bigger absolute value of a and (3 than the noisoLGM 
case. In this context, the information on the equation of 
state of the asymmetric nucleon-nucleon potential could 
be extracted by the temperature dependence of isoscaling 
parameters. 



in the lattice gas model. Here a' and (3' is the scaling 
parameters where we take Yi(N,Z) from 36 Ar instead 
of 36 Ca as shown above. Principally, the association of 
a number density with the exp(Ap/T) only is valid for 
a classical gas of free particles. The connection in this 
context of the LGM where the particles are interacting 
is not clear. However, as an attempt, we will still use the 
above equation to deduce the the relative free neutron 
density and proton density. In Fig. E{a) we show the p n 
(solid symbols) and p p (open symbols) as a function of 
N/Z of the systems at T = 5 MeV for the isoLGM case 
(squares) and the noisoLGM case (circles). In both cases, 
nearly linear relations of p n and p p have been observed 
with the increasing of N/Z. The Fig. OJb) shows the ra- 
tio p n / pp. The solid (dashed) line represents the deduced 
Pn/ P P from solid and open squares (circles) of Fig. [SJa) 
with isoLGM (noisoLGM) calculation and the dotted line 
is the initial value of p n j p p which is calculated assum- 
ing neutrons and protons homogeneously distributed in 
a volume proportional to the nucleon number. In isoLGM 
case, the values of p n /p P for neutron-rich nuclei ( 36 S and 

Si) is much larger than the initial value and while the 
values of p n /Pp for proton-rich nuclei ( 36 Ca) is much less 
than the initial value. While, in noisoLGM case, the val- 
ues of p n / ' Pp also increases with the initial N/Z value 
of sources, which is basically originated from the initial 
difference of isospin of hot emitters. The stronger en- 
hancement of Pn/Pp in the isoLGM case may indicate 
of a neutron enrichment while a proton depletion in the 
nuclear gas phase. In this context, it may be consistent 
with the isospin fractionation effect which is a signal ex- 
pected for the liquid gas phase transition in asymmetrical 
systems ISSSElllS- However, this interpretation 
is not unique, since the larger p n / p v can be also directly 
attributable to the interaction and thus probably not an 
increase in neutron density. 

In summary, the isoscaling is investigated using the 
fragment yield from the equilibrated thermal source with 
the same mass number but different N/Z which was pre- 
pared by the lattice gas model with the asymmetrical po- 
tential between like nucleons and unlike nucleons. The 
isotopic scaling and isobaric scaling are observed for light 
clusters and the isoscaling parameters a and (3 are ex- 
tracted from the ratios of R2i(N,Z) for fixed proton 
number or neutron number. It is found that a and -(3 
is almost the same and they drop with the temperature. 
However, the difference of the neutron (Ap n ) and pro- 
ton (A/ip) chemical potential does not change much with 
temperature and a slight and wide valley for (A/z ra — Ap p ) 
is observed around T = 5 MeV even though we have 
larger error bars for lower and higher temperatures, may 
indicate the onset of phase separation into isospin asym- 
metric dilute and isospin symmetric dense phase where 
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the liquid gas phase change occurs. The relative free neu- 
tron density or proton density is attempted to be deduced 
from the isoscaling parameter and they reveal a nearly 
linear relation to the N/Z of the initial system. However, 
the values for neutron-rich source are much larger than 
the initial value of N/Z as well as the values for neutron- 
poor source are much less than the initial value of N/Z. 
This may be from the isospin fractionation effect, or can 
be directly attributable to the interaction. 

In addition, in order to investigate the effect of the 
asymmetrical nucleon-nucleon potential in LGM, we also 
adopt the same potential between like nucleons and un- 
like nucleons in LGM. In this case, isoscaling still remains 
but the temperature dependence of isoscaling parameter 
(a and j3) vanishes and their absolute values decrease. 
The insensitivity of the isoscaling parameter to temper- 
ature means that the isospin partition between different 
fragments is almost the same regardless the excitation 
extent of system and the chemical phase separation is 
absent in this case. However, for the isoLGM, the asym- 



metrical nucleon-nucleon potential between like pair and 
unlike pair gives an additional contribution from isospin 
dependent equation of state to the isoscaling behavior ex- 
cept from the trivial contribution from the initial differ- 
ence of the isospin of the source, which results in a bigger 
absolute value of a and j3 than the noisoLGM case. In 
this context, information on the equation of state of the 
asymmetric nucleon-nucleon potential could be extracted 
by the temperature dependence of isoscaling parameters. 
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